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1. Introduction 

The current Tevatron collider and the upcoming Large Hadron Collider (LHC) experiments 
require a good understanding of the Standard Model signals to carry out a successful 
search for the Higgs particle and physics beyond the Standard Model. At these hadron 
colliders QCD plays an essential role and modeling all of the aspects of the events is crucial. 
From the lessons learned at the Tevatron, we need fixed order calculations matched with 
parton shower Monte Carlo's and hadronization models for a successful understanding of 
the observed events. 

Fixed order calculations are the first step in modeling the events. In this paper we 
describe an algorithm for the automated calculation of one- loop amplitudes. A successful 
implementation of numerical algorithms for evaluating fixed order amplitudes needs to take 
into account the so-called complexity of the algorithm. That is, how does the evaluation 
time grow with the number of external particles. An algorithm of polynomial complexity 
is highly desirable [1]. Another consideration is the suitability of the method within a 
numerical context. In particular algebraic methods can be successfully implemented in 
efficient and robust algorithms. This can lead to rather different methods from what one 
would develop and use in analytic calculations. 

Leading order parton level generators are well understood. These have been con- 
structed using algebraic manipulation programs to calculate the tree amplitudes directly 
from Feynman diagrams [2-6]. However, such a direct approach leads to algorithms of 
double factorial complexity. Techniques such as helicity amplitudes [7-10], color order- 
ing [11, 12] and recursion relations [13-17] have been developed and successfully used both 
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in analytic and numerical calculations of leading order amplitudes. These techniques all 
aim at further factorizing the calculation in smaller subsets. Of these the recursion relation 
technique stands out as it maximizes the factorizability of tree amplitudes and consequently 
can be implemented as a polynomial complexity algorithm of rank four [18] . Furthermore, 
recursion relations have a simple algebraic structure ideal for numerical evaluation of the 
tree amplitudes. Such a recursive algorithm was successfully implemented in tree amplitude 
generators [19, 20]. 

Compared to leading order generators the status of next-to-leading order generators is 
far less advanced. Calculations based on explicit one- loop Feynman diagrams using tensor 
form factor reductions have been used successfully for four, five and six point amplitudes. 
Examples of recent explicit calculations of six point one-loop amplitudes using such brute 
force methods are six photons [21, 22], six gluons [23] and e'^'e" four fermions [24, 25]. 
These direct approaches suffer from worse than factorial complexity. 

Alternatively, powerful analytic methods have been developed for one-loop amplitudes 
based on generalized unitarity methods [26-45]. These methods have been used in analytic 
calculations leading to compact expressions for processes such as for example vector boson 
plus four partons [46, 47]. The analytic methods developed are again based on the principle 
of factorization. By using generalized unitarity one can factorize the one-loop amplitude 
into tree amplitudes from which the coefficients of the master integrals can be constructed. 

The main problems of converting the on-shell analytic method into numerical algo- 
rithms are twofold. The first issue is resolving overlapping contributions between quadru- 
ple, triple and double cuts. Using the methods developed in [48, 49] one can disentangle 
the contributions in an algebraic manner. Secondly, the numerical implementation of the 
dimensional regularization within an on-shell method leads to complications. The on-shell 
cut lines give intermediate on-shell particles in non-integer number of dimensions. These 
lines are ill-defined and require special considerations. In particular from the aspect of a nu- 
merical implementation one needs a well-defined method for calculating the contributions 
originating from dimensional regularization. Several methods to calculate these contribu- 
tions applicable for numerical implementation have been developed [34, 35, 50-52]. These 
developments now allow the construction of numerical algorithms for evaluating one-loop 
amplitudes with a large number of external particles [53, 54]. 

In this paper we present an algorithm of rank nine polynomial complexity to calculate 
one-loop amplitudes with an arbitrary number of particles based on the purely algebraic 
methods developed in refs. [49, 51]. As a first step we implement the algorithm in the 
program Rocket^ using three and four gluon vertices to calculate all the helicity amplitudes 
of the pure gluonic scattering amplitudes. 

The outline of the paper is as follows. In section ^ we review the method and outline 
the construction of the required orthonormal basis vectors and D-dimensional polarization 
vectors. Section ^ contains our results: we study the accuracy of the method, the power-like 
growth of the computation time and give explicit single event results for various helicity 
configurations for up to twenty gluons. Finally, in section |^ we draw our conclusions and 

""^From the Italian Rucola: Recursive Unitarity Calculation of One-Loop Amplitudes. 
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give an outlook on our future plans. 



2. Construction of the Algorithm 

We implement the methods developed in refs. [49, 51] with some minor modifications into 
the Rocket program. These methods build upon the formalism of ref. [48] by removing 
the requirement of the four dimensional spinor language, thereby allowing for the extension 
of the method to D-dimensional cuts. 

To calculate the full one-loop AT-gluon amplitude, it is sufficient to be able to calculate 
the leading color ordered one- loop amplitude. From these color ordered amplitudes the 
full one- loop amplitude can be constructed [26, 55]. Eventually the one- loop amplitude 
has to be contracted with the tree amplitude, summed over the colors and integrated over 
phase space. The number of different orderings one needs to calculate can be drastically 
reduced by noting that the phase space integration will symmetrize the final state gluons. 
Alternatively, a more physical approach can be employed by not only randomly sampling 
over the helicities but also over the colors of the external gluons. In the following we will 
therefore focus on the leading color ordered amplitudes ylj^'(l, 2, . . . , N). 

We will use the master integral basis decomposition derived in ref. [51]. This decom- 
position in an overcomplete set of master integrals makes the loop momentum dependence 
on the dimensionality explicit: 

N / J 2 211213*415 21*2*3*426 

[n|*5] 

, V (/°'°^ T^^^ - rV(2,0) j{D+2) {D-A){D-2) (4,0) \ 
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[n|*2] ^ ^ 
where we introduced the short-hand notation [iili^j ] = 1 < ii < ^2 < • • • < in < and 

ID . = [J^ \ (2 2) 

withdj = di{l) = {l+Qi)^ = {l+Pi-\ hPi)^. In this basis the coefficients b\°fj, bf^fj, c\°fji^, 

4S2I3' ^^2^2324 ' 4lfXu' 4tSi3U^ and effXi^i^ are independent of the loop momentum 
dimensionality D. Because the coefficients d!-^fj^^-^, c-^fj^^ and 6^^^ are multiplied with 
a dimensional factor (D — 4) they cannot be determined using four dimensional cuts^. 
Therefore we need to extend the dimensionality of the cut line to higher dimensions. We 
choose the dimensionality to be an integer, resulting in a well defined on-shell particle after 
performing the cut [51]. 

^Note that the terms proportional to e^^f^^^i^^^ and d^^'il\^i^ will contribute only at 0{e). 
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By applying quintuple, quadruple, triple and double Ds-dimensional cuts (where Dg > 
D denotes the dimensionality of the spin-space) we can determine the coefficients of the 
parametric form of the one-loop amplitude. This requires the calculation of the factorized 
un-integrated one-loop amplitude 



Resj-^.-.j^^ (^y (0) 



X dij^j X A^jJ 



(0 



M 



{Ai,...,Am}=1 \k=l 



;{Afc) 



+1' 



_J^i^k+l)\ 
«fc+l ' 



(2.3) 



where M < 5 and the D-dimensional loop momentum / has to be chosen such that di-^ (I) = 
• • • = dij^,j{l) =0. As a result of the on-shell condition, the tree amplitudes ^'''1 have in 
addition to the external four dimensional gluons, two Ds-dimensional gluons with complex 
momenta. Dimensional regularization requires that Dg > D such that the ultra-violet 
poles are regulated. These higher dimensional gluons have [Dg — 2) polarization states. To 
calculate these tree amplitudes we use the standard Berends-Giele recursion relation [13] 
which is valid in arbitrary dimension and for complex momenta. 

The generic solution for the loop momentum in eq. (|2.3D is given by 



21 ■••«M «1---«M 



a 



M 



+ 




(2.4) 



for arbitrary values of the variables (Xi? 

The vector Vi^...^^ is defined in the space spanned by the denominator off'set momenta 
{gii, . . . ,qij^j}, while the orthonormal basis vectors {n^,j, ■ ■ ■ span the space orthog- 

onal to the space spanned by these momenta [49, 51]. Given the solution to the on-shell 
conditions in eq. ( p. 41 ), the loop momenta flowing into the tree amplitudes /j^, and 

in eq. ( p.3| ) are fixed by momentum conservation (see ref. [49]). 

Finally we need a stable and general method for constructing the orthonormal set of 
{D — M) basis vectors and the {Dg — 2) polarization vectors. To accomplish this we use 
the generalized Kronecker delta tensors [56] given by 



^Vl ^U2 ■ ■ ■ 



Jui Uu2 



... 5\ 



P2 
I'R 



(2.5) 



We use the notation 



= XlJ-llJ-2---lJ.R^ „ 



f2 



The [R — l)-particle Gram determinant is then given by 



(2.6) 



A(/ci, A;2, • • • , kn-i) = 5\ 



felfc2---fc_R-l 



(2.7) 



^The conformal transformation applied on the coefficients ai is only valid when V^^...^^^ 0. For the 
special case Vj^^...i^j = a similar solution is found. 
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Note that for R > D + 1 the generahzed Kronecker delta is zero. For the special case 
R = D we have the factorization of the Kronecker delta into a product of Levi-Civita 
tensors- A'^i'^^-'-md _ _/ii/i2 -7ii3p 

Given a set of momenta {qi, q2, ■ ■ ■ , qm} in a D-dimensional space-time we describe here 
how to construct the orthonormal set of basis vectors {ni, . . . ,nD-A4} such that Qi ■ rij = 
and rii ■ rij = 6ij. The set of momenta span the M-dimensional sub-space. The basis vector 
set {rii} spans the orthogonal {D — M)-dimensional space. For the A:-th basis vector we 
choose the arbitrary vector b^- The vector is then given by 

^/J. bf,_i---biqi---qM 

'^k ~ — ' \^-°) 

VA(6fc_i, . . . . .,qM)A{bk, ...,bi,qi,.. ..qu) 

with qi ■ riK = and nj ■ Uk = ^ik- Note that arbitrary vector ba-M in the basis vector 
^D-A/ drops out trivially because (Jjli-..^^ — ^Mi -'AtD ^ ■■^o. This means that in the 
construction of the basis we used (D — M — 1) arbitrary vectors. 

We can also use the above construction of the orthonormal basis vectors to define the 
Ds-dimensional polarization states of the cut gluon lines. Given a Dg-dimensional gluon 
with light-cone momentum p we want to construct a set of the {Dg — 2) polarization vectors 
{e^*^}. The polarization vectors have the property p ■ e*-*) = 0, e^*^ • e^^^ = —6ij and the spin 
sum is 

E = + ^^^^^^i^ -'Z^. (2.9) 

^ ^ ^ p - V [p- v)'^ 

where is the polarization gauge vector. Note that if is a light-cone vector the last 
term on the right-hand side is zero. The polarization vectors are easily constructed. We 
construct the {Dg — 2) orthonormal basis vectors with respect to the two four-vectors 
{p, f}. We now define the {Dg — 2) polarization vectors as e^P = irij . Then p ■ e^*) = 0, 
e(«) . g(i) = —fi - . = —Sij, and eq. ( ^ ) is satisfied. 

Once all coefficients in eq. (|2.1| ) have been determined using an appropriate set of cuts 
and loop momentum solution vectors, we can algebraically continue the dimensionality to 
the non-integer limit: D — > 4— 2e. This limit can be performed in different manners, leading 
to different schemes. In the four dimensional helicity scheme [57, 58] the continuation is 
defined as Dg — > 4, L) — > 4 — 2e while Dg > D. In the 't Hooft-Veltman scheme [59] the 
limit is defined as ^ 4 - 2e, D ^ 4 - 2e while Dg > D. 

Because we are interested in the next-to-leading order contributions we can neglect 
terms of order e in the continuation of the dimensionality. We then find for the color 
ordered one-loop amplitude 

^ j(0,0) .(4-2e) V;^ (0,0) .(4-2e) >p ,(0,0) j{A-2e) 

[n|«4] [iilia] [n|i2] 

j(4,0) (2,0) 2 

_ ^ "■iii2isi4 _|_ ^ ^1213 _ ^ Wi ~ g»2 j ^(2,0) _|_ ^^^^j ^2 10) 

[il\u] [«l|«3] [«l|«2] 
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The terms in the first hne give rise to the so-cahed cut-constructible part of the ampli- 
tude [60]. The terms in the second line can be identified with the rational part. In the 
approach used in this paper the division between these two contributions is not relevant. 

For the numerical evaluation of the bubble, triangle and box master integrals we use 
the package developed in ref. [61]. 



3. Results 

The algorithm is implemented in the Rocket Fortran 95 program. We first perform a 
series of checks by comparing to existing analytic results in the literature and by performing 
internal consistency checks. Next the accuracy of the results up to eleven external gluons 
is examined by comparing to the analytically known one-loop A^-gluon helicity amplitudes. 
We then show that the evaluation time as a function of the number of external gluons is 
consistent with a degree nine polynomial. Finally, we present explicit results for a few fixed 
phase space points for selected helicity configurations with up to twenty external gluons. 

In next-to-leading order calculations, the external gluons of the one- loop amplitude 
are identified with the momenta of well-separated jets. We therefore impose the following 
set of cuts on the generated phase space momenta 



hi < 3 , p±,^ > 0.01 V^, R^j = ^4 + > 0.4, (3.1) 

where r/j andp_L^j denote the rapidity and transverse momentum of particle i, (pij = {(pi — cpjl 
and r]ij = ~ Vjl denote the distance in azimuthal angle and rapidity between particle i 
and j and s = {Ei + £^2)^ is the center of mass collision energy squared. 

The unrenormalized one-loop A^-gluon results are given in the four dimensional helicity 
scheme (FDH). The conversion to the 't Hooft-Veltman (HV) scheme is given by 

^v,HV = ^v,FDH - -^Aree , (3.2) 

with 

{AttY r(i + 6)r^(i-6) 
= 16^2 r(i - 26) — • ^^-^^ 

We will give explicit results for double poles, single poles and the constant part after 
extraction of the factor cr- The pole structure of the A/^-gluon one- loop amplitude is given 
by [62-64] 



\i=i 



Av,poles = cr ( + 7 ( l^ln ( ) - ^ ) I ^tree, (3.4) 




where as usual the indices are modulo N. For the dimensional scale /i we use the center of 
mass collision energy, = s. 

3.1 Checks of the Results 

We have performed several checks on our numerical implementation: 

• we checked that the pole structure of the one-loop amplitudes agree with eq. (| 
up to twenty gluons; 
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• the number of loop momentum solutions to the on-shell conditions in eq. ( p.4| ) is 
infinite. This means we can solve for the coefficients in the parametric form of the 
integrand using many different sets of equations. We verified that the results are 
independent of the chosen loop momentum solutions; 

• we checked that the results are independent of all auxiliary vectors introduced in e.g. 
eq. ( |2.8| ) used to construct both the orthonormal basis and the polarization vectors; 

• we verified that our results are independent of the choice of dimensionality of the cut 
lines; 

• for six gluons we compared all one-loop helicity amplitudes with the numerical results 
of [23]; 

• for up to twenty gluons we compared with the analytically known one-loop helicity 
amplitudes where all gluons have positive helicities [33, 65]; 

• for up to twenty gluons we compared with the analytically known one-loop helicity 
amplitudes where all but one gluon have positive helicities [33, 65]; 

• for up to twenty gluons we compared with the analytically known one-loop helicity 
amplitudes where all but two adjacent gluons have positive helicities [26, 60, 66]. 

All checks performed were successful. 
3.2 Study of the Accuracy 

To study the numerical accuracy of the on-shell method implemented in Rocket we define 

I 1 VjUnit ^v,anly i 

where "unit" denotes the result obtained with the on-shell method and "anly" the analytical 
result for the constant parts of the one-loop helicity amplitudes (or in the case of = 6 
the numerical results of [23]). Similarly, where relevant, we denote by ^dp and esp the 
accuracy on the double and single poles, respectively. 

In fig. (H) we show the accuracy for the two adjacent minus helicity gluon MHV one- 
loop amplitudes, A^^{ !-■■■+)) for N ranging between six and eleven. The 100,000 

phase space points used for each multiplicity are generated uniformly in phase space using 
the Rambo algorithm [67]. We plot the accuracy for the double pole (X = DP [dp], solid, 
red), the single pole (X = SP[dp], green, dot-dashed) and the constant part (X = C[dp], 
blue, dotted). We first examine the six-gluon plot (fig. (|^), top left) and see that an 
excellent accuracy can be reached for all three contributions, the position of the peak 
being at eop = lO'^^.s^ ^ 10"^^-^, and ec = lO"^"-^, respectively. The tail of the 
distribution reaching to large values of e contains only a very few points for the single pole 
and the constant term. This lack of agreement is due to numerical instabilities. The well- 
known sources of instabilities are vanishing Gram determinants or other small intermediate 
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Figure 1: Accuracy on the double pole, single pole and constant part of the MHV amplitude with 
adjacent negative helicities for 6 up to 11 external gluons. Double ([dp]) and quadrupole ([qp]) 
precision results for 100,000 phase space points are shown. See text for more details. 



denominators. Several techniques have been developed in the past do deal with such 
exceptional points, such as developing systematic expansions [68-70] or interpolating across 
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the singular regions [71]. Similarly to what is done in [5, 50, 54], we adopt here a more 
brute force approach and recur to quadrupole precision. In (fig. (||), top left) we see three 
more curves: they correspond to the numerical accuracy on the same 100,000 phase space 
points when the one-loop amplitude is computed in quadrupole precision.^ One sees that 
the positions of the peaks move even more to the left, the peak of the double pole is now 
at edp = 10~^^'^ (magenta, dot-dashed line, labelled X = DP[qp]), of the single poles is at 
esp = 10~^^'^ (light blue, dot-dashed, labeled X = SP[qp]) and of the constant part is that 
ec = 10-13-2 (black, dot-dashed, labelled X = C[qp]). More importantly, out of 100,000 
phase space points samples, not a single point has an accuracy worse than lO"^. 

We can now examine what happens when the number of external gluons is increased. 
At double precision we can see from the position of the peaks and the width of the dis- 
tributions that the accuracy slowly worsens with increasing N. This is due to a slow 
accumulation of errors when more terms are added together and to the fact that there 
are potentially more instabilities. However, at quadrupole precision we see no appreciable 
worsening of the accuracy with increasing A'^. For = 11 the peak of the double pole is 
now at eDP = lO^^^'^, of the single poles is at esp = lO^^^'^ and of the constant part is that 
ec = lO-i^'*^. Again, out of 100,000 phase space points sampled, not a single point has an 
accuracy worse than lO"''. Up to A = 11 (and probably even for more gluons) quadrupole 
precision is sufficient to guarantee an accuracy needed for any next-to-leading order QCD 
correction. If higher precision is desired one can choose to evaluate the few phase space 
points which have insufficient precision using an arbitrary precision packages such as [72], 
at the cost of higher computation time. 

We note that while the plots here presented are for the MHV amplitudes, we performed 
a similar study for the finite amplitudes {A^^{+ ■ • • +), A^j!j{ — h • • • +)) and obtain very 
similar results. This indicates that the accuracy is essentially independent of the helicities 
of the external gluons. 

For the results shown in the above plots we choose to rerun all events in quadrupole 
precision to get an overall picture. However, in practice only a small fraction of phase 
space points require a quadrupole precision treatment (this fraction can be read off the 
plots in fig. (1^) and depends on A and on the target accuracy). Therefore one needs a 
systematic procedure to decide which events should be re-evaluated in quadrupole precision. 
One possible way is to verify the accuracy of the single poles results. The analytic single 
pole result is given in eq. (^^) for arbitrary number of gluons. Since two-point functions 
contain single poles, this checks the coefficients of the two-point master integrals as well as 
the coefficients of the higher point master integrals. In fig. (||) we investigate the correlation 
between the accuracy of the single pole contribution and the constant part. We plot the 
relative accuracy of the constant part log^gC^c) versus the accuracy on the single poles 
log]^o(^sp) in double precision (left) for A^ = 6 MHV amplitudes when sampling 1,000 
phase space points. The high correlation between the accuracy of the constant part and the 
single pole is evident. In fig. (^) (right) we show the improvement when running the same 
points in quadrupole precision (note the different scale on the y-axis). This leaves us with 

*Only the coefficients of the master integrals are computed in quadrupoie precision, master integrals are 
still calculated in double precision. 
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Figure 2: Relative accuracy of the constant part versus tiie accuracy on tlie single poles in double 
precision (left) and quadrupole precision (right) for A*" = 6 MHV amplitudes when sampling 1,000 
phase space points. 

a straightforward estimate of the accuracy of the one-loop evaluation. By comparing the 
numerical evaluated single pole result against the analytic single pole result we can decide to 
switch to quadruple precision to re-evaluate the full one-loop amplitude and get the required 
precision. Alternatively, one might choose to run in quadrupole precision whenever any 
small denominator (for instance in the construction of dual vectors) appears. We did not 
fully investigate yet which method is the most efficient in detecting potential instabilities. 
We will leave this to a future study, but we anticipate that identifying dangerous phase 
space points in order to increase the accuracy is not an issue. 

3.3 Time Dependence of the Algorithm 

Compared to traditional Feynman diagram approaches, the power of this method is that 
the time needed to compute one-loop amplitudes does not grow factorially with the number 
of external legs. It is indeed quite straightforward to estimate the scaling of time with the 
number of gluons N. Within the so-called constructive implementation of Berends-Giele 
recursion relations (or alternatively a recursive implementation "with memory") the time 
required to compute tree level ordered amplitudes grows as rtree,N oc N'^ [18]. Altogether 
the number of tree amplitudes that one needs to evaluate at one-loop is simply given by 

ntree = { {Ds, - 2f + {D,, - 2f ] (3.6) 



N\ (N\ N\ 

X 5C5 ,max ,max ,max o +2C2 ,max 



N 



where the first factor is due to the sum over polarization of the internal cut gluons in Dg^ 
and dimensions respectively needed to determine the dimensional dependence of the 
one-loop amplitude [51]. The constants Cm,max denote the number of times one needs to 
perform a multiple cut in order to fully constraint the system of equations determining 
the master integral coefficients. Explicitly we have cs^max = Ij C4,max = 5, C3^max = 10, 
and C2,max = 10 [51]. The integer number in front counts the number of tree amplitudes 
per multiple cut, finally the binomial coefficients corresponds to the number of possible 
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Figure 3: Time in seconds needed to compute tree (blue, dashed) and one-loop (red, solid) ordered 

amplitudes with gluons of alternating helicity signs, ^5^'(H 1 l----)i ^ function of the number 

of external gluons ranging between 4 to 20 using a single 2.33 GHz Xeon processor. 

cuts (for two point functions we subtract the vanishing contributions of the external self 
energy graphs). From this it follows that the time needed to evaluate a one- loop ordered 
amplitude will for large scale as 

Tone-loop,N ~ '^tree • ^tree.N OC . (3.7) 

In fig. dH) we plot the time needed to compute tree (blue, dashed) and one-loop (red, solid) 

ordered amplitudes with alternating helicity signs for the gluons, ^J^^(H 1 — . . .), as a 

function of the number of gluons in the range between four and twenty. Time estimates 
refer to using a 2.33 GHz Xeon processor. One can see that the times needed to compute 
tree and one-loop ordered amplitudes are consistent with a A^^ and growth respectively 
(we show a polynomial fit to those points as well).^ When running in quadrupole precision 
rather than in double precision the evaluation time grows with a factor of approximately 
thirty. We verified that the scaling with N is unchanged. 

Finally we remark that the time needed to compute "easier" or more "difficult" helicity 
amplitudes is approximately the same, i. e. we checked that the plot looks essentially 
identical for other helicity configurations. 

3.4 Results for fixed phase space points 

In this section we present sample results for one-loop helicity amplitudes at fixed phase 
space points for = 6, 7. Other results for N ranging from eight to twenty are given in 

^The evaluation times given in ref. [54] for 6, 7 and 8 gluon MHV amplitudes make use of analytic 
expressions for the required tree amplitudes. Nevertheless, the average time for the most complicated 6 
gluon helicity amplitudes quoted is 72ms per phase space point, while we have an average computation 
time of 90ms for any of the helicity amplitude. 
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Appendix |A[ Since the algorithm employed is independent of the chosen external helicity 
vectors we show, apart from a comparison with known amplitudes, the results for the "most 

difficult", alternating sign helicity configurations, ^Jv^(H 1 — • • •)• ^^^o in the tables 

means results smaller than 10^^^.) 

= 6 We choose in this case the same phase space point as in [23], i.e. the six momenta 
Pi are chosen as follows, p = {E,px,Py,Pz),^ 

pi = (-3.000000000000000, 1.837117307087384, -2.121320343559642, 1.060660171779821) 

P2 = (-3.000000000000000, -1.837117307087384, 2.121320343559642, -1.060660171779821) 

P3 = (2.000000000000000, 0.000000000000000, -2.000000000000000, 0.000000000000000) 

P4 = (0.857142857142857, 0.000000000000000, 0.315789473684211, 0.796850604480708) 

P5 = (1.000000000000000, 0.866025403784439, 0.184210526315789, 0.464829519280413) 

P6 = (2.142857142857143, -0.866025403784439, 1.500000000000000, -1.261680123761121) . (3.8) 

We obtain the same results as in [23], which we give in Table 1 for completeness. 



Helicity amplitude 




cr/e 


1 


(+ + + + ++)! 

|A^'™'(+ + + + ++)! 
+ + + ++)! 




1.060660419488780 ■ IQ-^'^ 




3.813284749527035 ■ lO-^* 


1.767767365814634 ■ lO^^^ 
0.529806483643855 
0.529806483661295 


I^tree + + + + 

|A^'™*(- + + + + + )! 

+ + + + + )! 


1.011255761241711 ■ 10~" 
2.377895374324842 ■ lO'^^ 


6.753625348984687 ■ 10-^° 
8.549005883762705 ■ lO'^^ 


3.963158957208070 ■ lO"!* 
3.25996704351899 
3.25996705427236 


lAl'^^i- - + + ++)\ 
|A^'"'^'*(- - + + ++)! 
|A^'"™(- - + + ++)! 


170.947689902659 
170.947689902659 


614.590878376396 
614.590878376397 


28.4912816504432 
1373.74753500854 
1373.74753500828 


I^r" (- + - + -+)! 
|A^'™'(- + - + -+)! 
+ - + -+)! 


18.8322923700467 
18.8322923700485 


67.7058293474830 
67.7058292869577 


3.13871539500808 
151.043950328960 
151.043950337947 


(+- + -+-)! 
|A^'™'(+ - + -+-)! 
- + -+-)! 


18.8322923700554 
18.8322923700485 


67.7058292857048 
67.7058292869577 


3.13871539500808 
153.780101529836 
153.780101415986 



Table 1: Results for tree level and one-loop virtual (unrenormalized) amplitudes in the FDH 
scheme for some helicity configurations for the case of six external gluons for the phase space point 
of eq. (3.8). Comparison with known results is also shown. 



N = 7 We randomly choose the following phase space point: 

pi = (-3.500000000000000, 3.500000000000000, 0.000000000000000, 0.000000000000000) 

P2 = (-3.500000000000000, -3.500000000000000, 0.000000000000000, 0.000000000000000) 

P3 = (1.721020317835363, 0.501455810979379, 0.991016354865482, -1.314663298501285) 

P4 = (2.156731348769508, -1.341260229883955, 1.169853444323893, 1.218176516478751) 

P5 = (1.037618172453970, -0.152260900991984, -0.998573413613423,0.237316723937150) 

P6 = (0.584066733036484, 0.353785124568635, 0.163503402186405, -0.435013415594666) 

P7 = (1.500563427904675, 0.638280195327925, -1.325799787762357, 0.294183473680051) 

The results are given in Table 2. 

®Note however that we use a different convention to label the momentum components here. 
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Hclicity amplitude 




cr/e 


1 


+ + + + + +)! 

+ + + + + +)! 

+ + + + + 







1.256534542409480 ■ 10"^° 
1.250170111559883 ■ lO'i^ 




0.310169532972026 
0.310169533483183 


\Atf'=<={- + + + + + +)] 

+ + + + + +)! 
lA^-'^^'yc- + + + + + +)! 


3.678212874319657 ■ 10"" 
2.713533399763100 ■ 10"" 


7.209572152581734 ■ 10"" 
8.924875144594874 ■ lO^^^ 




0.192052814810810 
0.192052814765395 


- + + + + +)! 

- + + + + +)! 

- + + + + +)! 


14.7462898421614 
14.7462898421614 


48.5008939631214 
48.5008939631213 


2.10661283459449 
87.3152155138790 
87.3152155138651 


+ - + - + -)! 

+ - + - + -)! 
+ - + - + -)! 


0.771305976661093 
0.771305976661095 


2.53684348996073 
2.53684348996075 


0.110186568094442 
5.93361050294547 
N.A. 


- + - + -+)! 

- + - + -+)! 

+ - + 


0.771305976661093 
0.771305976661095 


2.53684348996074 
2.53684348996075 


0.110186568094442 
6.04201240991614 
N.A. 



Table 2: Results for tree level and one-loop virtual (unrenormalized) amplitudes in the FDH 
scheme for some helicity configurations for the case of seven external gluons for the phase space 
point of eq. (|3.9| ). Comparison with analytical results, when available, is also shown. 



4. Conclusions and outlook 

In this paper we present the numerical implementation of an integer dimensional on-shell 
method for calculating one-loop amplitudes. The method used was developed in refs. [49, 
51] which we followed closely. The resulting program, Rocket, is a Fortran 95 code. The 
only limitation on the number of external particles are the available computer resources. 

As a first example and test of the implemented algorithm we calculated purely gluonic 
color ordered one- loop amplitudes. We explicitly study the properties and time-behavior 
of the algorithm up to twenty external gluons. The scaling of the computer time needed to 
evaluate color ordered one-loop amplitudes is consistent with the theoretical estimate of a 
rank nine polynomial. Comparisons to existing analytic results for special helicity config- 
uration shows a good numerical accuracy and only for a limited set of events quadrupole 
precision is required. For completeness and later reference we give numerical results for 
explicit events. 

Now that the algorithm has been validated up to twenty external particles we plan to 
include internal and external (massive) quarks and external vector bosons. This allows us 
to compute one-loop amplitudes to a large range of processes relevant for both the LHC 
and Tevatron experiments. 

We envision the collection of one-loop matrix elements to be integrated through match- 
ing into parton shower Monte Carlo's. The parton shower Monte Carlo will integrate the 
one-loop matrix elements with the radiative contributions. This will allow a seamless in- 
tegration of the available tools in a single framework and will provide a complete and 
advanced analysis tool for experimentalists. 
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A. More numerical results 

In this appendix we list some explicit results for high multiplicity events. 
= 8 We randomly choose the following phase space point: 

pi = (-4.000000000000000, 4.000000000000000, 0.000000000000000, 0.000000000000000) 

P2 = (-4.000000000000000, -4.000000000000000, 0.000000000000000, 0.000000000000000) 

P3 = (1.449692512284710, 0.958721567196264, 0.596442789329140, 0.909239977027169) 

P4 = (1.009340416556955, -0.511164560265637, 0.004406444973267, 0.870321464785572) 

P5 = (1.065731003301513, 0.640468049755497, -0.851783840257115, -0.006894789139810) 

P6 = (1.402207989626767, 0.338467194877458, -0.023915684844839, -1.360534911049079) 

P7 = (1.702524230799814, -0.823031265059939, 1.480876536483203, -0.167967189914708) 

PS = (1.370503847430242,-0.603460986503644,-1.206026245683657,-0.244164551709144). (A.l) 



The results are given in Table 3. 



Helicity amplitude 


cr/e2 


cr/e 


1 




A|''<=<=(+ + + + + + ++)! 
A^'™'(+ + + + + + ++)! 
A^-''"'^ (+ + + + + + ++)! 




3.853462894343397 ■ lO^^^ 




1.441159379540454 ■ 10"" 




0.196700600695691 
0.196700600738201 




A|''°<=(- + + + + + ++)! 
A^'™'(- + + + + + ++)! 
A^'''"'^(- + + + + + ++)! 




1.805821909003967 ■ lO"!* 


1.965638104048654 • lO^^" 
6.753606439965886 ■ 10"" 


2.257277386254959 ■ 10" 
0.528774716493063 
0.528774717652170 




^8''°°(- - + + + + ++) 

A^'™'(- - + + + + ++)! 

Al''"^'^(-- + + + + ++)\ 


34.6655135573561 
34.6655135573568 


129.645805347145 
129.645805291409 


4.33318919466960 
274.299773434926 
274.299773434900 




(-+-+-+-+)! 

A^'™'(- + - + - + -+)! 
Al'^'^'^i- + - + - + -+){ 


0.580921809110773 
0.580921809110846 


2.17259368023597 
2.17259368244769 


7.261522613885579 • 10"^ 
5.47630381976679 
N.A. 




Al"=<=i+- + - + -+-]\ 

A^'™'(+ - + - + -+-)! 
A^'""'^(+- + - + -+-)| 


0.580921809110862 
0.580921809110846 


2.17259368781042 
2.17259368244769 


7.261522613885579 ■ 10"^ 
4.92550054630729 
N.A. 



Table 3: Results for tree level and one-loop virtual (unrenormalized) amplitudes in the FDH 
scheme for some helicity configurations for the case of eight external gluons for the phase space 
point of eq. ( A.l ). Comparison with analytical results, when available, is also shown. 



= 9 We randomly choose the following phase space point: 

pi = (-4.500000000000000, 4.500000000000000, 0.000000000000000, 0.000000000000000) 

P2 = (-4.500000000000000, -4.500000000000000, 0.000000000000000, 0.000000000000000) 

P3 = (0.837513535184208, -0.699991554911619, -0.341245468610472, 0.308208168000110) 

P4 = (0.173971130750340, -0.127314067058472, -0.036582419646361, 0.112777698311325) 

P5 = (1.527678295320022, -0.329915424080283, -1.066446319295822, -1.042904135098816) 

P6 = (1.087863978393825, -0.427338530001958, -0.959598549771965, 0.282843489474564) 

P7 = (2.837576052736588,0.906765089674632,2.135787640996821, -1.633409342380760) 

P8 = (2.065394708000926, 0.474995872769939, 0.683997980573425, 1.890074332734359) 

P9 = (0.470002299614090, 0.202798613607760, -0.415912864245626, 0.082409788959218). (A.2) 



The results are given in Table 4. 



Helicity amplitude 




cr/e 


1 




A*'-<=<=(+ + + + + + + + +)! 
A^'™'(+ + + + + + + + +)! 
A^'""'^(+ + + + + + + + +)! 


4.860269836292316 ■ 10^^^ 
2.693624076029116 • 10"^^ 


1.845193695700690 ■ lO^o* 
1.176695244346755 ■ lO^^^ 


2.992915640032351 ■ lO"^"* 
5.66655561706295 
5.66655558047311 




Al"=<=i- + + + + + + + +)\ 

A^'™'(- + + + + + + + +)! 
A^^""'^(- + + + + + + + +)! 


3.938371378126140 ■ 10"" 
8.202679137223861 ■ 10"" 


2.340429860576292 ■ lO-o^ 
3.583296428617654 ■ lO'^^ 


9.114087930248735 ■ 10" ^"^ 
1.06208646061428 
1.06208646798175 




Al''=<=(- - + + + + + + +)] 

A^'™'(- - + + + + + + +)! 
A^'''"'^(- - + + + + + + +)! 


290.906701096922 
290.906701150957 


1270.81033486132 
1270.81033630185 


32.3229667945508 
3625.43061670521 
3625.43061670594 




+-+-+-+-)! 

A^'™*( + + + +)! 

A^'''"'^(- + - + - + - + -)! 


4.08169769666186 
4.08169769731155 


17.8306776720814 
17.8306776807844 


0.453521966367950 
57.1063950462874 
N.A. 




Al-^o" {+- + - + - + -+)] 

A^'""'*(+ - + - + - + -+)! 
A^'""'^(+ - + - + - + -+)! 


4.08169769662055 
4.08169769731155 


17.8306776454842 
17.8306776807844 


0.453521966367950 
55.0153807707576 
N.A. 



Table 4: Results for tree level and one-loop virtual (unrenormalized) amplitudes in the FDH 
scheme for some helicity configurations for the case of nine external gluons for the phase space 
point of eq. ( A.2 ). Comparison with analytical results, when available, is also shown. 



= 10 We randomly choose the following phase space point: 



pi 

P2 
P3 
P4 
P5 
P6 
P7 
P8 
P9 
PlO 



-5.000000000000000, 
-5.000000000000000, 
1.532520310665362, 1. 
0.258067002946760, - 
0.933615908667822, - 
1.406380992739491, 0. 
0.826140594064319, 0. 
1.003938368361987, - 
0.626103283169757, - 
3.413233539384504, - 



5.000000000000000, 0, 
-5.000000000000000, 
513070043279136, -0, 
0.183609750404875, - 
0.324162340665056, - 
945471671828880, -0, 
546878011685178, -0, 
0.166638867622269, 0, 
0.367431900674052, 0, 
1.963576867426942, 0. 



000000000000000, 
0.000000000000000, 
187195319404816, 
0.180504757414873, 
0.045792611756259, 
641345241905772, - 
615447980566360, - 
778081660907425, - 
140872898707371, - 
751331351433283, 2 



,000000000000000) 

0.000000000000000) 

155548896254710) 

0.017437606394811) 

-0.874334305926935) 

0.820162846752305) 

0.068238586686725) 

0.612137782060901) 

0.486984543874634) 

,688871562651979). 



(A.3) 



The results are given in Table 5. 



Helicity amplitude 




cr/e 


1 




+ + + + + + + ++)| 
^l6™'(+ + + + + + + + ++)! 
^I6""'''(+ + + + + + + + ++)! 


2.616999209810146 ■ 10"" 
7.645214091184737 ■ 10^^^ 


7.453142378465002 ■ 10"''^ 
3.853184186191476 ■ lO^^^ 


7.645214091184737 ■ lO^i* 
18.4349011284670 
18.4349011284671 




^!r(- ++ + +++ + ++)! 
^I6™'(- + + + + + + + ++)[ 
+ + + + + + + ++)! 


1.729567134060808 ■ 10"" 
3.138928592085274 ■ lO-^^ 


3.462486730362966 ■ lO"'"^ 
1.582018484813023 ■ IQ-^^ 


3.138928592085274 ■ 10"^^ 
14.1180690283674 
14.1180690283692 




AioH- - + + + + + + ++) 

^ir''(- -++++++++)! 
^ir"'(- -++++++++)! 


4899.72695665607 
4899.72695666341 


24694.6000400099 
24694.6000476827 


489.972695666341 
75844.9101458089 
75844.9101457814 




A-YoH- + - + - + - + -+) 
+ + + + +)| 
+ + + + +)| 


93.4611371998759 
93.4611372008902 


471.043678702711 
471.043677247939 


9.34611372008902 
1481.27447605664 
N.A. 




A*^f''(+ - + - + - + -+-) 
^l6™'(+ + + + + )l 

^ir"'(+ -+-+-+-+-)! 


93.4611371995618 
93.4611372008902 


471.043674005742 
471.043677247939 


9.34611372008902 
1503.97025803111 
N.A. 



Table 5: Results for tree level and one-loop virtual (unrenormalized) amplitudes in the FDH 
scheme for some helicity configurations for the case of ten external gluons for the phase space point 
of eq. (|A.3|). Comparison with analytical results, when available, is also shown. 
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= 15 We randomly choose the following phase space point: 



pi 




(-7.500000000000000, 


P2 




(-7.500000000000000, 


P3 


= 


(0.368648489648050, 0. 


P4 


= 


(0.985841964092509, - 


P5 


= 


(1.470453194926588, - 


P6 


= 


(2.467058579094687, - 


P7 


— 


(0.566021478235079, - 


P8 




(0.419832726637289, - 


P9 




(2.691168687878469, 1. 


PlO 




(1.028090983779864, - 


Pll 




(1.377779821947130, - 


Pl2 




(1.432526153404585, 0. 


Pl3 




(0.335532948820133, 0. 


Pl4 




(1.085581415795683, 0. 


Pl5 




(0.771463555739934, 0. 



7.500000000000000, 0.000000000000000, 0.000000000000000) 
-7.500000000000000, 0.000000000000000, 0.000000000000000) 
161818085189973, 0.125609635286264, -0.306494430207942) 
0.052394238926518, -0.664093578996812, 0.726717923425790) 
0.203016239158633, 0.901766792550452, -1.143605551298596) 
1.840106401193462, 0.715811527707121, 1.479189075734789) 
0.406406330753485, -0.393435666409983, -0.020556861225509) 
0.214182754609525, 0.074852807863799, -0.353245414886707) 
868400546247601, 1.850615607221259, -0.571568175905795) 
0.986442664896249, -0.193408556327968, 0.215627155388572) 
0.155359745837053, -1.074009172530291, -0.848908054184264) 
621168997409793, -0.290964068761809, 1.257624811911176) 
244811479043329, 0.138986808214636, 0.182571538348285) 
330868645896313, -0.756382142822373, -0.704910635118478) 
630840621587917, -0.435349992994295, 0.087558618018677). 



(A.4) 



The results are given in Table 6. 



Hclicity amplitude 






1 


\A\Ti+ + + + ---)\ 

Kr"(++++---)i 

|A--'^(+ + + + ...)! 














1.07572071884782 
1.07572071880769 


+ + + + 

\Air'\- +++■■■+ +)\ 
iAir'''(- +++••■++)! 














0.181194659968483 
0.181194659846677 


l^irC- - + + + ■•• + +)! 

i^ir'''(- -++■■■ ++)i 


111.867315217633 
111.867315217633 


586.858955605213 
586.858955605213 


7.45782101450887 
1810.13038312828 
1810.13038312852 


\Alf H- + -■■■ + -)\ 

\Air\- + ---- + -)\ 
iAir'"(- +-•■•+-)! 


8.776559143021942 • 10"^ 
8.776559143233895 • 10"^ 


0.460420629357800 
0.460420661976678 


5.851039428822597 ■ lO^^ 
1.52033417713680 
N.A. 


iA!r(+ -+•■■-+)! 
\Air\+-+----+)\ 
i^ir'^(+ -+■■•-+)! 


8.776559143021942 ■ 10"^ 
8.776559143233895 ■ 10"^ 


0.460420565320471 
0.460420661976678 


5.851039428822597 ■ 10"^ 
1.52960647292231 
N.A. 



Table 6: Results for tree level and one-loop virtual (unrenormalized) amplitudes in the FDH 
scheme for some helicity configurations for the case of fifteen external gluons for the phase space 



point of eq. (A^). Comparison with analytical results, when available, is also shown. The present 
results have been obtained by running in quadrupole precision. 
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= 20 We randomly choose the following phase space point: 



pi 






P2 







P3 




(0.325540699096246, - 


P4 


= 


(0.848759269032669, - 


P5 




(1.570982650317940, 0. 


P6 




(0.553167263468587, - 


P7 




(0.503893441998193, - 


P8 




(1.342531690799994, - 


P9 


= 


(2.116396457930369,0. 


PlO 




(0.602748352923314, - 


Pll 




(1.497270156443179, - 


Pl2 




(1.403440100824101, 0. 


Pl3 




(1.968885859795150, 1. 


Pl4 




(0.537434314204394, 0. 


Pl5 




(2.339779276321334, - 


Pl6 




(0.894504093025053, - 


Pl7 




(0.306394773317926, - 


P18 




(0.560235842911576, - 


P19 




(1.070093313805907, 0. 


P20 




(1.557942443784069, 0. 



10.000000000000000, 0.000000000000000, 0.000000000000000) 
-10.00000000000000, 0.000000000000000, 0.000000000000000) 
0.312416065575525, 0.033104753113790, 0.085305474968826) 
0.024590847182333, -0.848289178399890, 0.013894488596740) 
440773273547400, 0.720246141672893, -1.324745599853646) 
0.230303833105897, -0.480156399163221, -0.149679651832391) 
0.080227334603420, -0.333549326254518, 0.369075903611134) 
0.248744151369669, 0.500198592589739, 1.220786244980225) 
026327408340262, -0.849622753763036, -1.938190395961762) 
0.444717464828695, -0.325013925746484, 0.244740477851469) 
1.005250846935450, 0.200161613080160, -1.091432079774134) 
989809299131399, 0.670534066899850, -0.735054918411486) 
683278410651605, -0.987433952876609, -0.260882176167652) 
448194620535905, -0.219499970243363, 0.199441688897247) 
1.532082130972216, 0.290600646542316, 1.744374578491559) 
0.210636900721777, 0.805755963464840, 0.326384735908015) 
0.289150435585858, 0.066148943112216, -0.076773042418471) 
0.343584006920049, 0.286594014336811, 0.337161831792864) 
323329544655228, 1.003388464788536, -0.183764236276963) 
809991460939089, -0.533167693154030, 1.219356675598426) 



(A.5) 



The results are given in Table 7. 



Helicity amplitude 




cr/e 


1 




^2r(+ + + + ■■■)! 

^26"""(+ + + + ■■•)! 

^26'"'"(+ + + + ...)l 












1.78947750851720 
1.789477509283 




^*2r(-++ + - ■• + +)! 
Air\-+++---++)\ 

Air''i- + + + ■■■ + +){ 












0.303337144522210 
0.303337141901917 




^2r(- - + ++•■■ + +)! 
^2r"(- -++•■•++)! 
Air''i--++---++)\ 


334.192303003683 
334.192303003683 


1995.15970325579 
1995.15970325579 


16.7096151501841 
6882.49682704505 
6882.49682704481 




^2r(- +-■•• + -)! 
Air\- +-■■■ + -){ 

^26^"'"(- + -... + -)l 


4.194124200605681 ■ lO""* 
4.194124200039273 ■ lO""* 


2.503931965487835 ■ 10"^ 
2.503931873702081 ■ 10"^ 


2.0970621000196 ■ 10"^ 
8.456871985787404 ■ lO'^ 
N.A. 




A%--{+ - + ...- +)\ 
Air'i+- + ---~+)\ 
^26'"'"(+- + ...-+)l 


4.194124200605681 ■ lO"" 
4.194124200039273 ■ 10""' 


2.503931902332899 ■ IQ-^ 
2.503931873702081 ■ 10"^ 


2.0970621000196 • 10"^ 
9.203156962017870 ■ lO'^ 
N.A. 



Table 7: Results for tree level and one-loop virtual (unrenormalized) amplitudes in the FDH 
scheme for some helicity configurations for the case of twenty external giuons for the phase space 



point of cq. (A_^). Comparison with analytical results, when available, is also shown. The present 
results have been obtained by running in quadrupole precision. 
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